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Abstract 

Theoretical work predicts that the strong dependence of T c on pure shear 
strain within the a-b plane of optimally doped YB^CusO/.^ results in het- 
erogenous columnar pins of vortex lines about dislocation lines and about 
nano-columns inclusions aligned in parallel to the c axis. The critical cur- 
rent of a rigid vortex lattice driven by the Lorentz force in the presence of 
such clusters of pin/antipin lines is computed using two-dimensional (2D) 
collective pinning theory and by numerical simulation of the corresponding 
2D vortex dynamics. Both theory and computer calculation find that the an- 
tipin component of the heterogenous columnar pins contributes substantially 
to the net in-field critical current. 
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1. Introduction 

Optimally doped films of the high-T c superconductor YBa 2 Cu 3 07_5 (YBCO) 
grown by pulsed laser deposition (PLD) show high critical currents that reach 
a substantial fraction of the theoretical maximum depairing current in ap- 
plied magnetic field oriented in parallel to the c axis[l|. Evidence exists that 
this is primarily due to edge dislocations that thread the film along the c 
axis in PLD-YBCO[2j. The critical current of YBCO films in c-axis aligned 
applied magnetic field can also be enhanced substantially by the introduction 
of BaZr03 (BZO) or of BaSn03 (BSO) material that naturally form nano- 
rods along the c axis during the film growth process^ [ij] (5)]. Such increases 
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in the in-field performance of the critical current in YBCO films by columnar 
microstructures are very important to the development of second-generation 
high-T c superconducting wire technology. Achieving a fundamental under- 
standing of how they occur is therefore desirable. 

In this paper, we shall compute the in-field critical due to natural or to 
artificial columnar microstructures in YBCO thin films by employing two- 
dimensional (2D) collective pinning theory^ [3] and by performing direct 
computer simulations of the corresponding 2D vortex dynamics 



Pinning of vortex lines will be assumed to be due to variations in T c about 



columnar defects that in turn are a result of the strain field about them 15 



Heterogeneous columnar defects are thereby predicted that are composed 



of component lines of pins and of antipins of equal strength [16|[17j. These 
form a dipolar "clover- leaf" pattern in the case of an edge dislocation line 
defect (see fig. [I]), while they form a quadropolar pattern in the case of a 
nano-column inclusion (see fig. [2]). Both 2D collective pinning theory 0] [3] 
and direct computer simulation of 2D vortex dynamics obtain the following 
central result: that the critical current due to an arrangnement of weak het- 
erogeneous columnar pins is equal to that expected from an arrangement 
of conventional columnar pins[l8|] obtained by unbinding the component 
pin/antipin lines from the heterogeneous cluster, and by then converting 
the isolated antipin lines that result into pin lines of equal strength. Figure 
[3] depicts the two pinning landscapes in question. The antipin components 
of heterogeneous columnar defects therefore increase the critical current by 
a factor of two compared to the contribution due to the true pin components 
alone in the collective pinning regime. 

2. Strain/pinning landscape about material line defects in YBCO 

It is well known that the critical temperature in optimally doped YBCO 



couples most strongly to pure shear strain in the a-b plane [151] . In particular, 
the variation in the critical temperature depends linearly on the symmetric 
strain tensor e Qj/ g = {d a up + d/3U a )/2 following 

ff c (r) = ^^(aT c /M^(r) (1) 



with strain derivatives dT c /de aa = 230 K and dT c /debb = — 220 K. Below, we 
shall show how the shear strain about a line of edge dislocations or about a 
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nano-column inclusion in a YBCO thin film combined with this experimental 
fact results in a heterogeneous columnar pinning line for a vortex line in the 
mixed state. 

Consider a line of edge dislocations that threads a YBCO film in parallel 
to the c axis. Gurevich and Pashitskii calculated the variation in the critical 
temperature expected from a linear dependence on strain ([I]) and obtained 



the result 16 
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with constants! 1 

1 / dT c dT c \ 

BT BT 

T ± = ^-^ = -450K. (4) 

Above, a denotes the Poisson ratio and b denotes the magnitude of the 
Burgers vector. It can be aligned along either the a (+) or the b (— ) axes. 
We now note that the contribution of the vortex core to the vortex line 
tension is approximated by the fundamental energy scale per unit length e — 
($ /4vtAl) 2 , where Xl denotes the London penetration depth, and where <3>o 
denotes the flux quantum in a type-II superconductor. The potential-energy 
landscape experienced by a vortex line then has a contribution fei(r) = 
{deo/dTc) ■ 5T c {r) due to the variation in the critical temperature about the 
dislocation line ([2]). It is dominated by the shear component (Tj_) in the case 
of optimally doped YBCO, which is shown in fig. [TJ The potential energy 
landscape is essentially dipolar in such case, with two anti-pinning lobes (+) 
and two pinning lobes (-). 

An analogous heterogeneous pinning potential for a vortex line about a 
nano-column inclusion has been obtained more recently by one of the au- 



thors in collaboration with Barnes and Varanasi 17]. The lattice mismatch 
between the inclusion and the YBCO matrix results in pure shear strain 
about the inclusion. These authors then ultimately find a variation of the 
critical temperature in the YBCO matrix about the nano-column given by 

(in) 

tTJ.r)=T ±7 L „„;■ (^) 2 CO3 20. (5) 

^out C | + C\ ' V r 7 
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Here Ar = r in — r out is the mismatch between the radius of the inclusion, r in , 
and the radius of the corresponding void in the YBCO matrix, r out , while 
(m,out) an( ^ c (m,out) ^ e 2D compression and shear moduli, respectively. 
The potential-energy landscape experienced by a vortex line again then has a 
contribution 5ei(r) = (deo/dT c ) -ST c (r) due to such a variation in the critical 
temperature about the nano-column inclusion, but with d-wave symmetry. 
Figure [2] depicts the constant-energy contours. Notice that the potential 
energy is characterized by diametrically opposite lobes of pins and of antipin 
oriented with respect to each other by 90 degrees. 

Last, the force fi(r) = — Vei(r) that a vortex core experiences due to 
changes in T c about both a line of edge dislocations ([2]) and about a nanocol- 
umn inclusion §5§ is long-range. Suppose now that a random field of parallel 
lines of heterogeneous pins thread a YBCO film along c axis at transverse 
locations {R«}. Statistical cancellations in the net force experienced by a 
vortex core, f(r) = ^jfi(r — Rj), cut off the long-range force. The effec- 
tive finite range r p of each heterogeneous pin is determined by matching the 
fluctuation in the net force over the system with the long-range force due 
to a single material line defect: f 2 = |fi(r p )| 2 . It is given by the expres- 
sion r p = (2/7™^) 1 / 6 r^ in the case of strain-induced pinning of vortex lines 
about nano-column inclusions [17|, where denotes their density. 

3. Collective pinning of Abrikosov vortex lattice by heterogeneous 
columnar pins 

Above, we have shown how the strain field about edge dislocations and 
nano-column inclusions that thread a YBCO film along the c axis result in 
heterogeneous pinning lines for vortex lines that are present in the mixed 
phase. As shown by figs. (pQ) and ([2]), the heterogeneous columnar pins are 
clusters with equal numbers of pins and antipin, each of equal strength. [See 
also Eqs. ([2]) and (JSJ).] Consider now a YBCO film threaded by a field of such 
heterogeneous pinning lines in applied magnetic field aligned along the c axis 
as well. Next, suppose that a rigid vortex lattice appears as a result. Among 
the fraction of vortex lines that experience forces from the heterogeneous 
columnar pins, half will lie in valleys of the vortex-core potential while the 
other half will lie on the hills of the vortex-core potential. Because the 
magnitude of the force experienced by the vortex lines that lie in the valleys 
and on the hills of the vortex-core potential are equal, the critical current 
should be the same as that due to a field of columnar pins which results 
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from (?) unbinding the cluster of pins and antipins from the heterogeneous 
pins, and (ii) converting each of the antipins that are now isolated into pins 
of equal strength and range. Figure [3] depicts the two columnar pinning 
arrangments in question for the case of nano-column inclusions. Below, we 
shall see that this prediction is indeed borne out by 2D collective pinning 
theory. 

Previous work by one of the authors and MaleyQ implies that many het- 
erogeneous columnar defects collectively pin a rigid Abrikosov vortex lattice. 
Theoretical arguments and numerical Monte Carlo simulations indicate, in 
particular, that a hexatic Bose glass state can exist at low temperature in 
such Here the vortex lattice is threaded by isolated lines of 



edge dislocations parallel to the lines of relatively weak pins/anti-pins. (See 
fig. HI) The critical current j c is then limited by plastic creep of the vortex 
lattice associated with glide by such edge dislocations [3] . Application of 2D 
collective pinning theoryQ yields a transverse Larkin scale R c that is set by 
the density of Larkin domains 

R~c 2 ~ n p (f p /c m b A ) 2 . (6) 

Here n p denotes the density of vortex lines pinned by the material line defects 
and f p denotes the maximum force per unit length at a pin/antipin that 
makes up the heterogeneous material line defect. Also, c 66 denotes the shear 
modulus of the vortex lattice, while b A is the magnitude of Burgers vectors 
that glide. It can be shown that R c is of order the separation between isolated 
lines of dislocations in the hexatic Bose glass|21]. A statistical summation 
of the pinning forces inside of a corresponding Larkin volume, in which the 
vortex lattice is perfectly rigid, then yields a net pinning force density 

j c B/c ~ n p f p 2 /c 66 b A . (7) 

The above is valid if many vortex lines experience pinning forces inside of 
a Larkin volume: n p R 2 ^> 1. Last, study of equations (jSJ) and (jSJ) shows 
that the potential energy about a heterogeneous pinning line experienced by 
a vortex line has zero angle average. This implies that the occupation of a 
heterogeneous pin by a vortex line is purely random. The density of vortex 
lines that they collectively pin is then equal to n p = (<7pnB)n0, where de- 
notes the density of heterogeous columnar defects, where cr p = n(r p — r„ ut ) is 
the effective crossectional area of a heterogeneous pinning/antipinning line, 
and where % denotes the density of vortex lines. Observe now that both f p 
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and the product n ( f ) a p remain constant after unbinding the pins and anti-pins 
that make up heterogeneous columnar defects, and after then converting the 
isolated antipins that result into isolated true pins! We conclude that 2D 
collective pinning theory[3]j^] implies that the critical current due to com- 
parable arrangements of conventional columnar pins and of heterogeneous 
columnar pins like those shown in fig. [3] are the same. 



4. Computer simulation of 2D vortex dynamics 

We shall now compute the in-field critical current due to columnar pins 
instead by simulating over-damped vortex dynamics with the Langevin equa- 
tion. Our aim is to once again compare the critical-state performance of 
type-II superconductors threaded by conventional columnar pins[18( with 
that of type-II superconductors threaded by heterogeneous (d-wave) pinning 
lines that can result from the strain field about nano-column inclusions in 



YBCO|17l|. The Abrikosov vortex lattice in a field of parallel columnar pins is 



in a Bose glass state characterized by an infinite macroscopic tilt modulus 22 
The longitudinal Larkin scale L c is then infinite (l9|. which means that vor- 
tex lines may be considered as rigid rods at long longituginal length scales. 
In a film geometry that is perpendicular the applied magnetic field, their 
dynamics can therefore be described by the over-damped 2D Langevin equa- 
tions for the transverse motion of a given rigid vortex line at location r 
3 00 13 : 



VP — ^ F vv + F vp + F Lo rentz + F Br own(^)- (8) 
r' 

It is governed by the Bardeen-Stephen viscosity rj for a vortex line. Each 
rigid vortex line experiences forces F vv = —VU VV and F vp = —VU vp due to 
interactions with other vortices at locations r' and due to the field of material 
line defects, respectively. The former is taken to be long range: 

U (r r') = 4n£ ° d V exp[q-(r-r')] 

NxN y ^4-2 cos(g,a) - 2 cos(q y a) ' 1 ; 

Note that periodic boundary conditions are imposed over a rectangle of di- 
mensions N x a x N y a. The lattice constant a of the grid is understood to be 
of order the diameter of a vortex core. Above, d denotes the thickness of the 
film. Values for both potential energies U vv and U vp that span a N x x N y 



6 



grid are stored in look-up tables. The corresponding forces F vv and F vp are 
computed from these potentials by taking finite differences over the grid. 
They are also stored in look-up tables. Last, a force at an arbitrary location 
is determined by linear interpolation from those determined off the grid. 

In the special case of a featureless (S-wave) line pin that lies at a grid 
point (ma,na), the above linear interpolation scheme yields a pinning force 
of the form 

fx P) ( x , V) = Kd ( x -ma)-g(y- no) , 

fy p) ( x , y) = g( x ~ ma ) ■ fw(v - na )> 

where is the pinning force along a principal axis of the grid that contains 
the pinning center. Both it and its corresponding ID potential energy are 
shown in fig. [5j The scale /g that appears there denotes the magnitude of 
the maximum pinning force, which is related to the depth Uq of the line pin 
potential by /q P ^ = Uq /a. The transverse interpolation of forces is set by the 
function g, which is also shown in fig. It is given explicitly by g(y) = 1 — 
\y/a\ if \y\ < a, and by g(y) = otherwise. Notice now that d x fy P ^ ^ d y f x p \ 
which means that the above force field for a line pin is not conservative. On 
the other hand, leaving the ID function g free while imposing a conservative 
force implies that it must be given by g(y) = — «^(y)//o instead. As 
shown by fig. 0, both functions have qualitatively similar shapes. They 
both also have the same area, f_°°dyg(y) = a. We believe, therefore, that 
the model for a linear pin displayed by the equations above is qualitatively 
correct. 

The rigid vortex lines also experience a uniform Lorentz force FLorcntz that 
push them across the width of the wire, in addition to a random Brownian 
force F Br own(£) that brings the vortex lattice into thermal equilibrium. We 
numerically evolved the above Langevin equation (jSJ) in order to describe the 
over-damped dynamics of 896 rigid vortex lines. The second-order Runge- 
Kutta method was employed. Values for the interaction forces and for the 
pinning forces to be stored in look-up tables were determined by taking finite 
differences over a 224 x 224 square grid that spans the system, and which has 
a lattice constant a that is understood to be of order the superconducting 
coherence length. The computer simulations were conducted in two stages. 
We first annealed the vortex lattice down to low temperature in the presence 
of a landscape of weak columnar pins. Figure [3] depicts the two pinning (anti- 
pinning) landscapes that were studied. All (component) line pins had equal 
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strengths, each with a maximum pinning/antipinning energy = /g a 
equal to 1% of the dominant interaction scale 2eqcI. The initial tempera- 
ture was taken to be very many times the Halperin-Nelson- Young melting 
temperature of the 2D vortex lattice, Tm D \ and it was successively halved 
until a very small fraction of Tm D ^ was reached. Two million time steps were 
taken at each temperature, with a time interval At equal to a few percent 
of the natural time scale rja^/ 'Aneod that is set by the dominant interaction 
energy among vortex lines, Eq. 0. Here a v = ($ /-B) 1//2 is the average sep- 
aration between vortex lines. The annealed vortex lattice was found to be in 
a hexatic vortex glass state characterized by long-range orientational order 



and unbound dislocations 19] 20] 23]. Only a fraction of the vortices in the 
lattice experience pinning forces. In the second stage, we drove the annealed 
vortex lattice at zero temperature with a uniform Lorentz force in the pres- 
ence of the same landscape of columnar pins. The driving force increased in 
300 steps from zero to about twice the critical force (current). Each step of 
constant Lorentz force lasted 200, 000 time iterations, At, and the Lorentz 
force was increased at each step by a part in 100, 000 of the interaction scale 
2e d/a v . 

Figure El shows the current- voltage properties that result from driving the 
vortex lattice at zero-temperature through the pinning landscapes that are 
shown in fig. [3j The landscape in the left panel is the field of conventional 
(S-wave) columnar pins that is shown in the left panel of fig. [3j All of the 
pins in this case are featureless and of equal strength, with a number that 
matches the number of vortex lines (896). As depicted by fig. each pin 
lies at a grid point and inherits the minimum range of the grid. Notice the 
critical current that separates a pinned vortex lattice at low currents from a 
vortex lattice that experiences flux creep and flux flow at high currents. It is 
interesting to note that although the vortex lattice was threaded by isolated 
dislocation defects throughout the simulation, these became well seperated 
precisely at the critical state. (See figs. HJandEl) The dislocation defects then 
gradually began to pair up and annihilate with increasing driving force in the 
flux-creep/flux-flow regime. The right panel of fig. [6] displays the current- 
voltage property for a driven vortex lattice that experiences the arrangement 
of quadrapolar columnar pins that is shown in the right panel of fig. [3j The 
latter aims to model the strain-induced pinning potential shown in fig. [2] that 
is predicted to exists about nano-column inclusions in YBCO films 17| . This 
pinning landscape contains half the number (2 x 224) of true line pins, by 
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comparison with the previous case, with an equal number (2 x 224) of antipin 
lines. Each component line of pins/antipins has the same range and strength 
as in the previous case, on the other hand. Figure [6] indicates that the critical 
current in the present case of heterogeneous columnar pins is comparable to 
that obtained in the previous case of featureless columnar pins, which had 
twice the number of true pins! This result is nevertheless anticipated by 2D 
collective pinning theory, Eq. (j7j). It is due to the fact that the pin and the 
antipin components of the heterogeneous pinning line contribute equally to 
the critical current of a rigid vortex lattice. Last, two well-spaced dislocation 
defects with equal and opposite Burgers vectors threaded the vortex lattice 
throughout the second simulation. 



5. Conclusions 

In summary, strain about columnar microstructure in YBCO films gen- 
erates heterogeneous pin/antipin lines that can multiply the critical current 
expected from the true pin line components alone by a factor of two. This 
result is obtained for weak correlated pinning inside of the collective pinning 
regime, however. It is unlikely to persist in the limit of strong correlated het- 
erogeous pins, where all vortex lines that experience pinning forces then lie 
along a true pin component of the composite line defect. Antipin line com- 
ponents of the heterogenous pin have no effect on the critical current in that 
limit. A key assumption of the collective pinning analysis is that the critical 
state be found in a hexatic Bose glass state characterized by isolated disloca- 
tion defects that thread the vortex lattice 0]. Figure 0] demonstrates that this 
assumption is indeed confirmed by the molecular dynamics simulations that 
we conducted across the critical state for 896 rigid vortex lines. Previous 2D 
simulations find evidence for the break-up of the vortex lattice into grains at 
much higher numbers of rigid vortex lines when the latter experience short- 



range interactions |13j] [141] . Thermodynamic Monte Carlo simulations on over 
2000 vortices that experience the same long-range logarithmic interaction 
studied here, Eq. find evidence for an equilibrium hexatic vortex glass 
phase Ji| , on the other hand. It therefore remains to be seen if the hexatic 
glass phase obtained in the present molecular-dynamics simulations persists 
in the thermodynamic limit. 

JPR thanks Paul Barnes and Hsin-Ju Wu for discussions. This work was 
supported in part by the US Air Force Office of Scientific Research under 
grant no. FA9550-06-1-0479. 
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HETEROGENEOUS PINNING ABOUT EDGE DISLOCATION 




Figure 1: Shown is the contour of a potential-energy landscape for a vortex core about 
an edge dislocation aligned parallel to the c axes of optimally doped YBCO at two fixed 
energies ±e\. 
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HETEROGENEOUS PINNING ABOUT NANO-COLUMN INCLUSION 
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Figure 2: Contours for the potential energy of a vortex core about a nano-column inclusion 
in YBCO are depicted for two fixed energies ±£\. 
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896 PINS 224 *(2 PINS + 2ANTIPINS) 




Figure 3: Shown above are windows of equal area into an arrangement of conventional 
columnar pins (+) and into an arrangmcnt of quadrapolar columnar pins (+x, x+) that 
can result from shear strain about nano-column inclusions. 
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Figure 4: Shown is the critical-state configuration of 896 rigid vortex lines in the presence 
of an equal number of conventional columnar pins. (See the left panel in fig. [6l) Pairs of 
red and green dots mark dislocation defects in the triangular vortex lattice. 
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Figure 5: The solid lines above depict the pinning force and the corresponding potential 
energy along principal axes (x = or y = 0) due to a featureless point defect that lies on 
the square lattice grid at point (0, 0). The force at a general point (x, y) is determined by 
linear interpolation over the grid. 
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Figure 6: The left panel above shows the current- voltage curve obtained from driving 896 
rigid vortex lines in the presence of an equal number of conventional (S-wave) columnar 
pins via the Langevin equation ((8]). All forces and velocities are given in units of 2e d/a„ 
and of 2eod/r]a v , respectively. The right panel above depicts the current voltage curve 
obtained from driving 896 rigid vortex lines in the presence of 224 quadrapolar columnar 
pins. Each component pin/antipin line has a strength and range equal to that of the 
conventional S-wave columnar pins. (See fig. E3) 
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